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SUMMARY 


The present report consists of three AIAA papers concerning research 
carried out under NASA Grant NAG T-203 during the period 1 July 1982 to 30 
June 1982. Technical monitors on this research work were Dr. Douglas Price 
and Dr. Christopher Gracey of NASA-Langley's Theoretical Mechanics Branch. 

The papers are as follows: 

On-Board Rear-Optimal Climb-Dash Energy Management, A. R. Weston, 

E. M. Cliff and H. J. Kelley, presented at the America! Control Conference 
San Francisco, California, June 22-24, 1983. 

Optimal Symmetric Flight with an Intermediate Vehicle Model, P. K. A. 

Menon, H. J. Kelley :*nd E. M. Cliff, for presentation at the AIAA 
Guidance and Control Conference, Gatlinburg, Tennessee, August 15-17, 1983. 

Energy State Revisited, H. J. Kelley, E. M. Cliff and A. R. Weston, 
for presentation at the AIAA Atmospheric Flight Mechanics Conference, 
Gatlinburg, Tennessee, August 15-17, 1983. 


OF POOR QUAL8TV 




ON-BOARO fjEAR-0'»TIHAL CUHS-BASH ENERGV MANAGEMENT 
A, R, Weston^ 

E. H. C11ff2 
H. J. Kelley^ 

Aerospace and Ocean Engineering Department 
Virginia Polytechnic Institute and State University 
Blacksburgi Virginia 24061 


N84 16X16 


The subject of this paper is the study of optimal 
and near-optimal trajectories of high-performance 
fighter aircraft in symmetric flight. On-board, real- 
time, near-optimal guidance is considered for the 
climb-dash mission, using some of the boundary-layer 
structure and hierarchical ideas from singular per- 
turbations, In the case of symmetric flight this re- 
sembles neighboring-optimal guidance using energy-to- 
go as the running variable. However, extension to 
3-D flight is proposed, using families of nominal 
paths with heading-to-go as the additional running 
variable. Some computational results ?re presented 
for the symmetric case. 

Introduction 

On-board, real-time guidance is constrained by 
the limited computational resources available, par- 
ticularly on fighter aircraft, whore space and weight 
are at a premium. As a result the algorithms employed 
must be simple to implement and have small storage re- 
quirements. The objective of the presently described 
effort is to investigate practical algorithms for a 
variety of missions in 3-D flight. In the present 
paper anapproach is developed for the Intercept mission 
in synmetrlc flight based on a concept, sketched in 
Ref, 1, in which extensive numerical computation is 
required on the ground prior to the mission, but the 
on-board execution is simple. The scheme takes ad- 
vantage of the boundary-layer structure common in 
singular perturbations, studied in Ref. (2), arising 
with the multiple time scales appropriate to aircraft 
dynamics. Energy modelling of aircraft, as first 
examined in Refs. (3-5) and extensively developed in 
Ref. (6,7) is used as the starting point for the 
analysis. In the symmetric case, a nominal path is 
generated which fairs into the dash or cruise state. 
Feedback coefficients are found as functions of the 
energy-to-go, (dash energy less actual energy), along 
the nominal path, These serve to generate transitions 
towards the nominal path, closed-loop, and to counter 
disturbances. In this situation the guidance method 
is similar to the neighboring-optimal guidance methods 
of Refs. (8-161, However there are two significant 
differences! in the present work the gain Indexing 
is done in terms of the current energy; this avoids 
problems encountered in estimating the index time. 
Further extension to 3-D flight is considered here 
where families of reference paths would replace a 
single trajectory, with heading-to-go as the addition- 
al running variable. 
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Problem Formulation 


The overall objective is to develop an on-board, 
real-time near-optimal flight-control system in 3-D 
for a variety of missions and for arbitrary initial 
conditions. The example developed here is for the 
climb-dash Intercept mission in 2-D. The equations of 
motion for a point-mass model of an aircraft can be 
written: 


E=V(nT-D)/W (1) 
li “ V sin y (2) 
y = (L cos W cosY)/mV (3) 
X « V cos Y (4) 


The following assumptions are embodied: fixed mass, 
thrust along the path, flight over a flat earth, and 
no winds aloft. 


Aerodynamic Modelling 

The aircraft which is used as an example to per- 
form numerical calculations is a high-performance 
interceptor. The drag is modelled as a parabolic 
function of the control: 


The thrust is a function of Mach number and altitude; 
it is stored using the spline-lattice technique of Ref. 
17. The flight envelope is shown in Fig, 1. 


Nomenclature 


Reduced-Order Hodellin q 


Cjj Drag Coefficient 

Cjjp Zero-Lift Drag Coefficient 

^ Research Associate 
2 

Professor 


Order reduction, based on observed or assumed time- 
scale separations, is attractive in solving flight 
vehicle control problems. This is not only due to the 
smaller number of states and unknown initial conditions 
but nxire significantly to improvement of the condition- 
ing of the system of differential equations by the 
confinement of the more unstable dynamics to 
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boundary-1 ayot* corrections. It has been approclated 
since the work of R«f. 3 that h and y can be changed 
tnuch faster than E» which In turn Is a 'fast' variable 
In coniparlson to the range. This leads to the refer- 
nwlatlng of the equations of motion, following Ref. f, 
with the Interpoletlon parameters t' and |2 


fi * V sin Y (5) 

Y • (L - M Cos Y)/mV (6) 

c’ t • V(nT - D)/W (7) 

X » V cos Y (1) 


The Introduction of three separate time scales In the 
system must conform to the requirement of the Tihonov 
theory {Ref, 18ll that the ratio (c2/c') 0 as c' 0, 

as discussed In Ref. 6, Two resulting possibilities 
are now discussed. 

R ectnine ar- Hotlon Mod el 

The simplest model Is found when cl and are 
taken as 0. The result of these assumptions on the 
lillfferentlal equations are noted: 



« 0 E ' 0 nT»D (10) 

These equations embody the assumptions that the 
states h,Y. E can be varied Instantly In a control - 
like fashion, while n and Cjj, are chosen to satisfy (9) 
and (10), The control -like states h, £ are chosen to 
minimize the Hamiltonian, kx As Xx Is constant 
the m1n-H operation picks out the dash point on the 
envelope. This Is the zeroth-order "outer" solution 
of singular-perturbation theory, which the solutions 
from other time scales fair Into, 

Energy-State Models 

If cl “ 1, k 2 0 the next level of complexity, 
known aS energy modelling, Is found, y and h are 
still assumed to be *contro1-11ke‘ , and 'fast', while 
E Is now a 'slow' variable. Cj Is still determined by 
eq. (9) j h and n are chosen to minimize H^x^k + xee, 
where E Is governed by eq, (3) and Xe by 


n and h are found 1n terms of the ratio of X£ to X^ 
and their signs: the ratio determines the relative 
Importance of range rate and energy rate, while the 
signs determine the .s.’nsa of the optimization. For 
example If X£ Is negative, then n " 1 minimizes H‘, If 
Xe Is positive then n = 0. If Xx “ 0 h will maximize 
the excess power, leading to an 'energy-climb' (Refs. 

6, 19, 20) i the resulting path Is shown on the h-V 
plane In Fig. 2. This schedule shows multiple jumps 
In altitude which arise from realistic variations on 
the thrust data. This contrasts with other examples 
where the altitude jumps are typically due to the 
transonic drag-rise (Refs, 19, 20). 

In the Intercept problem x» f 0. The "energy climb' 
path traced out for Xx = 0 Is sfiown In Fig, 2, where 
the path does not fair Into the "outer" solution 
operation at the dash point. The analysis and com- 
putations are more complex: the value of xe Is now 
pivotal and the Initial value of (Xe/Xx) will de- 
termine the ensuing path. Finding the solution which 


fairs into the dash-point is therefore a two-polnt- 
boundary- value problem In one dimension, solved using 
the usual 1-D search techniques. This 'range-optlmal- 
energy-climb' 1s shown in Fig, 3, with the energy 
climb for comparison. 

The disadvantages of this approach are that It re- 
sults In a y approximation of 0 (whereas In the climb- 
dash problem the actual angle for a typical high per- 
foriimnce fighter may be 1r excess of 45°). and that 
Instantaneous jumps In the altitude appear at the end- 
points, and also are possible during the trajectory. 

Singular-Perturbation Pro c edure 

By the use of singular-perturbation theory, 
boundary- layer type correctUns can be used to over- 
con« the weaknesses of the energy model, I.e. Initial, 
Internal, and final jumps In altitude, While the 
altitude discontinuities are eliminated by expansion 
to the zeroth order, nonzero y values are obtained In 
the procedure of Ref, 6 by continuing the expansion to 
the first ordKi-. This Is a nontrivial problem In the 
case where the altitude transitions occur at the be- 
ginning or the end of a trajectory, and Is even more 
complex In the case of an Internal altitude Jump. As 
a result, the zeroth-order-corrected energy model loses 
Its attraction when realistic yS are required for on- 
board use as commands, 

On-Board Guidance 

An alternative to using order reduction, suggested 
In Ref. (l), Wi.lch is simple enough to lend Itself to 
on-board implementation Is now developed, for the case 
of symmetric flight, The scheme has roots in the 
hlearchical structure of solutions of the energy model. 
In which specific energy Is a relatively 'slow' vari- 
able and Its value determines the contrul-llke 'fast' 
variables, h and y. This suggests that trajectoiMes 
of the point-mass model funnel rapidly Into the vi- 
cinity of a single path, which fairs Into the dash- 
point. The idea pursued 1n this and the companion 
paper of Ref, 1 Is to determine this 'skeletal path', 
for the point-mass model, for as wide a range of 
energies as possible. This Is the nominal, or re- 
ference, trajectory and the altitude and path-angle 
histories are recorded as functions of the energy-to- 
go. The next step Is to generate a neighboring- 
optimal feedback guidance law which will control the 
aircraft so as to follow a neighbor of the nominal 
optimal path. Linear-feedback coefficients generate 
transients which bring the aircraft to the vicinity of 
the reference trajectory. The guidance law Is a linear 
feedback control based on the difference between the 
nominal and actual altitude and path-angle values. The 
coefficients, which correspond to minimizing the second 
variation, as 1n Refs, (8&16), are found by perturbing 
the altitude and path angle separately from their nomi- 
nal values along the reference trajectory, The optimal- 
control problem Isre^solved and the partial derivative 
of the control with respect to the states (at fixed 
energy) Is found In difference-quotient approximation. 
The Ct, commands which are sent to the autopilot are 
taken from the nominal path, with linear corrections 
for variation with altitude and path-angle from their 
nominal values: 

^ W' (h-h*(£)) + (E) (y-y*(E) 

On-board use for full-throttle climb-dash requires 
only the storage of the states h (£) y*(E), the con- 

3C| 

trol Cj*(E), and the feedback coefficients (E), 
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(E), as functions of the energy or energy-to-go, 


Optimal Solutions U sing T he Pol nt-Hass. M o del , 

A requirement of the proposed ides is a largo 
number of optimal -control solutiof s to the point-nv»is- 
modclled problem. This can be done by the use of 
direct methods, such as gradient methods, where the 
control history is parameterited by sectlonally-linear 
or spline approximation and the terminal conditions 
are met by either penalty or projectiyr. techniques, 
Alternatively the question can be posed as a two-point 
boundary value problem with split conditions and 
boundary conditions, which can be solved by indirect 
methods such as multiple-shooting (Refs. 21, 22). To 
solve the problem of time-optiital control, one forms 
the variational Hamiltonian: 

H - AgE + + a^y + A^j! 


and applies the Maximum Principle (Ref, 23, 24)). Tho 
resulting Euler differential equations are: 





iil 

ax 


The lift and the th ottle setting must be chosen (with 
the present sign convention) to minimize the Hamilton- 
ian, which requircjs that: 


optimal path using point-mass model dynamics, over the 
widest possible energy rang?. In the climb-dash 
problem being studied, the highest energy corresponds 
to that of tho high-speed point on the aircraft en- 
velope, the dash ‘outer solution'. The lowest energy 
corresponds to the trajectory which just kisses the 
terrain limit, i.e, below this energy, optimal so- 
lutions which start at zero attitude would dive below 
the terra'll limit if it were absent, 'his lower 
energy is found by examining the initial load factor 
of a family starting from level flight at the terrain 
limit altitude; when the initial load factor is 1» 
this lower energy is determined. This is shown in 
Fig, (4), where the initial load factor is plotted for 
several different initial energies. Further, as the 
climb is the boundary layer of the dash "outer" so- 
lution, it must fair into the dash state asymptotically, 
If a finite time is used, this is theoretically im- 
possible and a final control transient is inevitable 
to meet the bouridary conditions; however such tran- 
sients can be minimized by picking a climb trajectory 
which is long enough in duration. A measure of smooth- 
ness with which the final approach is made is the final 
load factor: as the time of flight is increased it is 
found that tho final load factor approaches a value of 
1, exponentially, see Fig, 5. 

A trajectory of sufficient duration that the final 
load factor is 1.001 is taken as the nominal path. This 
resembles an energy-climb schedule in that the energy 
of the aircraft determines the optimal altitude and 
path-angle; in fact the altitude chosen, for the Range- 
Optimal-Energy Climb is close to «,hat from the Euler 
Solution at the same energy; Fig. J compares the two 
solutions in the V-h plane. It is important to note 
that 28-dec1mal-d1git precision was required for '•jc- 
cessful convergence of the multiple-shooting itf native 
process for the "long" 2-0 point-mass Euler Solutions 
under discussion. 


and 
n • 1 
(Aj < 0) 


Method of Solution 

Euler solutions were found in the present effort 
by the method of multiple shooting, using the algo- 
rithm and computer program of Refs. (14, 25) kindly 
supplied by DFVLR, Oberpfaffenhofen, FR6, In this 
method, the interval of integration is broken up into 
many subintervals. This is preferable to 'simple 
shooting', as optimization problems of lifting flight 
are ill-conditioneo, the state-Euler system being in- 
herently unstable. Partitioning the time interval has 
the effect of suppressing numerical-error growth. 

This need arises in the calculation of the feedback 
gains, found by the difference of the control at the 
beginning of two optimal solutions. As a result, to 
find the gains to 4 figures the control must be known 
to at least 8 figures. The multiple-shooting method 
has greater accuracy than the other methods available, 
and although it is sometimes difficult to generate the 
initial reference trajectory, the subsequent calcu- 
lation of the feedback gains is rsJatively easy, as 
the method has good convergence properties in the 
vicinity of a solution. 

Optimal Reference Trajectory 
The first objective is to generate a reference 


Feedback Coe f*i .:-jnts 

The feedback copf ''i'rieiiti for the neighboring- 
optimal guidance schv. 5-‘e nv^dod as functions of the 
energy-to-go along the ,&ferencu path. As discussed 
earlier, these are found in difference-quotient ap- 
proximation as the average of the forwards and back- 
wards differences. The gains are found to be very 
large near the final state and the feedback system 
correspondingly 'nervous'. In the last foot of 
energy-to-go, both of the feedback gains start to grow 
exponentially. This sensitivity of neighboring - 
optimal-guidance schemes close to the terminal state 
has been noted in the literature (Refs. 8-16), 

Simul ation 

Following the satisfettory spHning of the nominal 
states, controls and feediiack coefficients as functions 
of the energy-to-go, the guidance scheme was tested by 
running a simulation of i,he point-mass-model , using the 
feedback law, and comparing the resulting trajectory 
with an Euler solution which started from the same 
initial conditions. With zero disturbance the auto- 
pilot was able to follow the nominal path more than 
satisfactorily, over the entire range of energies, de- 
spite the Inevitable errors which arise in the spline 
representations and numerical integration. Tests were 
performed with the initial altitude disturbed from that 
of the nominal path at different energies up to 15000 
feet above and below the nominal path. The resulting 
trajectories are shown in Figs. 7-12. These show that 
the feedback law follows the optimal solution closely, 
even when the initial disturbance is far outside of 
the range of linearity of the feedback gains. 
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Jnwlementatlon 7. 

The analysts so far has not taken into account the 
variation In aircraft welghtt 3-D flight, winds aloft, 
and nonstandard atmospheric properties. The problem B. 

must be re-solved with variable weights before actual 
on-board use. Also, winds aloft and temporaturo vari- 
ations could be accounted for by precomputing the 
variation In tho stored variables using first-order 9, 

approKimatlons. A method for enhanced aeroth-order 
approximation via special choice of state variables 1s 
described In Ref. 26 and appears compatible with the 
present approach, 10. 


Conclusions 

The numerical results b'tar out the following con- 
clusions: first, that all trajectories which fair in- 11. 
to the high-speed point consist of a rapid transition 
onto a reference or skeletal path. If they do not 
originate on It. Secondly, the ll.neav- feedback scheme 
proposed Is able to control the aircraft so that It 
closely follows the appropriate neighbor of the nominal 12. 
path for large perturbations of Initial conditions, A 
3-D extension of the computational scheme Is of In- 
terest in which there are two dominant states, I.e. 
headlng-to-QO In addition to energy-to-go. As a re- 13. 
suit, families of optimal paths which fair Into the 
dash-point will be needed, and the feedback coefficients 
will be functions of two variables (represented via a 
spline lattice) Instead of one, 14. 
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Fig, 1 Level -Flight Envelope 
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Fig. 6 Range-Energy CUmb and 

Corresponding Point-Mass Solution 
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Fig, 10 Altitude vs. Energy, 

Initiel Pointi5,000 ft. above 
NOffiinal 




Fig. 9 Altitude vs. Energy, 

Initial PointlO.OOO ft. above 
Nominal 


Fig. 12 Altitude vs. Energy 

Initial Point 10,000 ft, 
below Nominal 
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Abstract 

Optimal flight In the vertical plane with a 
vehicle model intermediate In complexity between 
point-mass and energy models Is studied, Flight- 
path angle takes on the role of a control variable. 
Range-open problems feature subarcs of vertical 
flight and singular subarcs as previously studied. 

The class of altitude-speed-ranqe-time optimi- 
zation problem', with fuel expenditure unspecified 
is Investigated and some Interesting phenomena un- 
covered. The maximum-! Ift-to-drag glide appears as 
part of the family, final -time-open, with appropria- 
te Initial and terminal transient maneuvers. A 
family of climb-range paths appears for thrust excee- 
ding level-flight drag, some members exhibiting os- 
cillations. Oscillatory paths generally fall the 
Jacobi test for durations exceeding a period and 
furnish a minimum only for short-duration problems, 

Minimizing paths of long duration follow a 
certain corridor in the V-h chart, The features of 
the family sharpen for the special case of thrust 
and drag Independent of altitude, and considerable 
analytical attention Is accorded to this for the 
Insight It provides to the more general model. 

The problem of "steepest climb" Is found to be 
m-posed with the vehicle model under consideration, 
stralght-vertkally-upward maneuver sequences being 
furnished by a family of paths alternating between 
upward and downward vertical flight and Including 
a limiting "chattering" member. 


Introduction 

There has been interest from the beginning of 
optimal-flight studies in approximations Featuring 
simplified vehicle models. Representation of drag 
as the drag for level flight leads to an Inter- 
mediate vehicle model In which path angle v takes on 
the role of a control variable and the order of the 
system 1s reduced by one, An additional order- 
reduction leads to an "energy-state" model with al- 
titude or speed as a control variable (Refs, 1, 2 
and 3). This 1s reviewed In a companion paper 
(Ref. 4). The present paper examines optimal sym- 
metric flight with the intermediate vehicle model, 

The analysis Is based in part upon an explora- 
tion of Euler solutions for the path-angle-as- 
control model carried out In Ref. 5. The present 
analysis examines higher-order optimality conditions 
and "chattering-control" phenomena. The weaknesses 
of the model will be seen as more extensive than 
previously noted. 

♦Research supported by NASA Langley Research Center 
under Grant NAG 1-203 

-Aerospace Engineering Graduate Student, Student 
Member AIAA 

-^Professor of Aerospace Engineering, Fellow AIAA 
--•^Professor of Aerospace Engineering, Member AIAA 


Intermediate Vehicle Model ; 

The point-mass dynamical model of aircraft 
flight incorporating the assumption of thrust-along- 
the-path Is given by 


V “ gC*^'^^'^ - sinY] 

0) 

fi “ Vslnv 

(2) 

A ■ VCOSY 

(3) 

w « Q 

(4) 

' « i (^ - cosv) 

(5) 


Here V is airspeed, h altitude, x down range, W 
weight fuel consumed, y flight-path angle, T 
thrust, D drag, g the acceleration due to gravity, 

L lift and Q the fuel-coniumptlon rate. 

» 

The sweeping asfumptlon that drag can be ap- 
proximated by Its level-flight value Is next 
Invoked. This permits the deletion of equation (5) 
and the elevation or path angle y to control status. 
Lift coefficient, Cl, or angle of attack a, pre- 
viously a control variable, 1s correspondingly 
assumed to be such as to satisfy (5), There 1s 
obviously trouble ahaad with this modelling should 
'r turnout to be large In optimized maneuvering or, 
worse yet, should y exhibit jump behavior. 

The optimal -control problem to be treated, 
then, Is the minimization of a function of the 
final values of the state variables and final time. 

The variational Hamiltonian function Is 
H • Sinv}* X^VSInv + X||VCos-, 

( 6 ) 

and the Euler-Lagrange equations are 

- ‘^V I § (7) 


Xy I yT-D) - I 


^x * ° 

\r\j ^ 0 

w 

and 


Xy gCosY + Xj^VCoSY - Xj^VSIpy “ 0 


(S) 

(9) 

( 10 ) 
( 11 ) 
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In the following, the t1 \e derivatives of equation 
(n) will be used to eliminate the time varying 
costates In favour of the control » and 
derivatives, Note that this Is somewhat formal 
since ■ may not exist. Using equations (7) - (11) 
one may now proceed to eT^fjInate those cpstates 
which are variable In the Hamiltonian. Using (11) 




( 12 ) 


and hence 


" i ( ^ ) + V ( L - X vSec^Y ) 

V - h^ X g " (13) 

substituting for v from (12) In (8) and using 
equation (7) 

X. + V( Xu- XJany ) H 19. « 0 (14) 

Using (13) In (7) and using equation (1), one ob- 
tains a second expression for 1}, 


x^, + x 


+ X. 


h I 


(7-0) + 


3 I 

w 


T-0) 






(15) 


Equations (14) ^e.d (15) may now be used to obtain 
an expression for Xj^ In terms of X^^ and xj;.. 


¥ 


5 V(T-D) - a 3 V(T-D) 

^ 7 W 


m 


- X, 


Tanv ^ 3V(T-0) ■ a 3 V(T-D) ',+ g 

tin 7 W vCosY 


- ySec^Y 


+ X\ 


r 




s 0 


(16) 


The expressions (12) and (16) may be used for vjllml- 
natlng and Xj, 1n the Hamiltonian with the 
following result; 


-Cosy 


■ X. 


Note that 


{If ’ ?■ tvj ^ ^ 


E = Constant 


Where E ■ h + V. 

2g 


, the Specific energy. 


In order to investigate the implications of this 
complicated sxpresslon, consider first the case of 
free final value of range x and fuel W. If the 
final values of these variables are left open, then 
the natural boundary conditions Jx“0 and “*^"0 apply 
and the optimization problem 1$ a trade-of“ between 
final values of time t, altitude h and airspeed V, 
the maximum or minimum value of one of these varia- 
bles or some function of these variables being 
sought without regard to range or fuel consumption. 

In equation (17), If the transversal Ity condition 
for minimum time, H ■ -1 Is Imposed, the well-known 
energy-climb schedule Is obtained. 

One notes that, In this case, equation (17) can 
be satisfied either by cosy « 0, vertical flight, 
or by vanishing of the bracketed expression, viz. , 
the partial derivative of specific excess power 
V(T-O) with respect to altitude with specific energy 
held constant. Thus the solution of this, or any, 
h-v-t optimum problem 1s made up of vertical climbs, 
vertical dives and "energy climbs" pieced together 
In the proper order. Similar considerations apply 
if fuel expenditure rather than time Is to be mini- 
mized. In this case H ■ 0, ‘x ’ ^ and * 1, and 
equation (17) yields the minimum-fuel-climb path 
with fixed throttle In the V-h plane. 

If range Is to be maximized or minimized with 
final time and fuel unspecified, then x^ » +1 and 
H » xm n 0, and a first-order differential equation 
for path inclination emerges as follows. 


’ 3J_ 

\ aF TW 


(7-0) - (T-D)V- r 
OO! 


:osy 


0 


(18) 


If one chooses x« » -1 and a fixed value of H 
(to be determined), with ■ 0, expression (17) Is 
the Euler equation for maximizing range to climb 
with fixed final time, With H ® 0, "• -1 and a 

fixed value of x.y, similarly, the maximum range to 
climb trajectory with fixed final value of fuel Is 
obtained. It may be noted that the maximum range 
to cUmb problem Is 111-rposed In that the range to 
climb for thrust greater than drag without time or 
fuel constraints does not have a maximum or even 
an upper-bound. Further, fixed-throttle range-fuel 
trajectories are not of any significant interest 
In practical situations. Hence, attention will be 
focused on the problem of maximizing the range-to- 
c11mb with a specified final time (fixed H, X;^ =■ -1) 

The system (1) - (3) and (18) generates a tra- 
jectory family for the range problem. The possibi- 
lity of obtaining analytical solution of the system 
for the case of thrust and drag as arbitrary func- 
tions of altitude and air speed is remote. However, 
using the assumption of constant-density atmosphere, 
wherein the thrust and drag depend on airspeed only, 
one can obtain an analytical solution to this 
system (Ref. 5). The expression (18) can be re- 
written as 


^ “ Tm {k ■ fl7 } 


(T-D) 


(19) 


Time may be eliminated in favor of airspeed V as an 
independent variable. 
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In the following, several tran'j formations of 
Independent variable are carried out without 
attention to monotonlclty requirements, the thought 
being to fit the solution segments ootalned Into 
families in due course. Tha temptation of range as 
Independent variable will oe avoided, hcn'ever, In 
anticipation of purely vertical motion segments, 

In the Interest of brevity we designate p » (T-0)/W 


With altitude-dependence suppressed, the path angle 
Y Is determined as the solution of the first-order 
differential equation 



Slny 


du 


( 21 ) 


Further simplification Is obtained oy another change 
of Independent variable, this time from V to p 



{ Slnv " u 


I 

w 


(22) 


If the roles of independent and dependent variables 
are now regarded as reversed, this equation takes 
the form 


2 

du + 1 u - n Slnv » 0 

dv Co7y Co'sy 


(23) 


which Is the form of the Bernoulli differential 
equation 


gbi + fi (y) + fp(Y) » 0 

cJy 


(24) 


with i = 1. According to Kamke, (Ref. 6), this 
equation has the solution 


L.£(v)y^dv 

where 

E(,) - ( 26 ) 


with Identification of f-j and fg as 

fl(v)* 1 

Cosy 

f,(v) . 

The solution (25) becomes as follows 
InCosY 

j = Cosy 
= SinY + C Cosy 



(27) 

(28) 


(29) 

(30) 


Before expressing this relationship in the form 
V = >(ii), we relate the Integration constant C to 
equilibrium values of u and y corresponding to un- 
acpelerated flight. Such values may be designated 
with a superscribed bar: 

u = SinY (31) 

C = CotY (32) 

The solution may then be expressed as: 


OF POOR QUALITY 


«,lnY Slnv + Cosy Cosy * g, 

ti 

or as 

Y » Y + COS^ j £ j 
Here u Is the value of w In unaccelerated flight and 

Y » Sln^u (35) 

In Fig. 1, the solution (34) Is Illustrated for 
various values of u. The range of angle y has been 
restricted to + 180® In this plot. 

With this solution at hand, ihe state histories 
can be generated. If the thrust Is taken as zero, 
the state»Euler system produces the flattest-gllde 
trajectory, flown with maximum llft-to-drag ratio, 
along with a family of transients to and from this 
point (Fig, 2). When a positive margin of thrust 
over drag exists, a family of oscillatory solutions 
is generated for various values of C as shown in 
Fig. 3. It may be noted In Fig. 3 that the Inner- 
most point corresponding to u * >2 in V-y space 
corresponds to flight at (T-O)fpax' while along the 
outermost closed path, the flight path angle v 
switches between t 90®. 


(33) 

(34) 


With the availability of the Euler solution 
(30) to the maximum-range problem with altitude 
dependence suppressed, one may proceed to obtain a 
similar solution to the more general Euler equation 
(17) using variation of parameters (Ref. 7). 
Equation (17) may be written as 


9 



CosS H 1 4Vg) 


+ (S ^ 3 


''x V^u 




m 


(36) 


As In equation (22), the Independent variable’ 
is changed from time to airspeed resulting In 


^ { SinY - W ) ’ Cosy du - Cos y H 
dV u 3V ■ 




Cos^Y 0 ^ d 


1 ^ 
X V-y 

V 


,(Vy) 




(37) 


Rearranging, one obtains 


dy( Siny - y 

av 


- Cosi3y = - CosS H _Ul,(Vy) 

y W ■ ^ 


^X ^ ^ (38) 


+ cos^Y ji d 


Equation (30) Is the analytical solution to 'the 
differential equation (38) with H and \(|, both zero. 
The expression (30) may be differentiated with 
respect to airspeed to obtain 


- 1 8y = ( Cosy - C Siny ) ^ Cosy (39) 

TW 3v 3v 

y 


Note here tnat C Is no longer a constant, but a 
function of the independent variable V. Substituting 
for w in (39) from (30) 


B3 


1 3u 

7 w 


- jicosY - C stnv ) ^ ^ 


[ sinv + t '^osy"1 


(40) 


Using equation (40) In (38) 

COS^Y ^ ‘ " Cos^Y H 1 d (Vu) 

SInY ■*■ C Cosy 37 ^ 

Since II * 1 from (30), 

STHT+l-Coi? 

Cos'^Y tiC » Cos y 

37 


(41) 




X ^'V/ 
tW 




du 


1 

V p 


SI] 


(42) 


The quantities within thei Ibrackets can be Identif- 
ied as ' ^ 


[rir] 


-d 

37 

and 

-d 

37 


From Which 


= 1 + _L_ dp 

vtj2 ^ 



1 . ^ + 1 ^ + 4 - 

7737 37^ 


n* 




VS 


dC = 

37 




Equation (45) is readily Integrated to yield 


C = 


tL 1 ' jS 4_ ^ ^1 

r 7T t 7T 


(43) 

(44) 

(45) 

(46) 


1 = Siny 

7 


where C, is an arbitrary constant. Hence for the 
time-range-fuel problem, the solution with altitude 
dependence suppressed is 

w L - ^ X * y»»'' (47) 

\ / 

To express the above result in the form v = 
•,(,.), we need to relate the integration constant 
Cl to equilibrium values of u and y corresponding 
to unaccelerated flight. Unlike the situation in 
the simpler problem, the Interpretation of equation 
(47) is not straightforward. 

From a practical viewpoint the time-range 
problem is of main interest since minimum-fuel 
problems with fixed throttle are rare. Fuel -range 
problem will not be discussed further in the present 
paper and in subsequent development the fuel mul- 
tiplier will be taken as zero. 

Investigation of equilibrium points with 
■iM = 0 results in a plot of the values of H/)v vs 
airspeed as shown in Fig. 4 for a parabolic (T-D) 
distribution illustrated in Fig. 5. In Fig. 4 
three separate regimes can be identified. H/Xx 
values to the left of the (T-O)niax velocity are 
positive while those between the (T-D)(nax 
the V(T-D)niax Point have a negative sign. All H/'x 
values to the right of the speed for V(T-D)niax 
positive. Any of these values may be used to 


evaluate the arbitrary constant ci as follows. As 
in (31) 

ORIGINAL PAGii m ^^8) 
OF POOR QUALIW 


■ Sin? 
V » V 


Equilibrium value of H, 


Coty - 1 L« 1 _ + C, 
^X V P 


or 


C, « Coty - IL J 
Xx V M 

using (51) in (47) 


(49) 

(50) 

(51) 


i * Siny 
u 


putting 4 » |_1_ - + Cotyjand using 

a well known Trignometric identity, 

Y * Tan^ rn + Cos’ f 1 1 

LJ 


(53) 


Equation (53) is the Euler solution to the 
time-range problem with altitude dependence of u 
suppressed. In Figs. 6, 7 and 8, the analytical 
solution evaluated for representative H/Xx values 
from each of the three regimes is shown. Fig. 6 
and 7 indicate oscillatory solutions and are in the 
neighborhood of a stable equilibrium point. The 
similarity of these figures to Fig. 3 is striking. 
The solutions in Fig. 8 are non-oscillatory and bear 
some resemblance to. Fig. 2. 

Summarizing, one notes that the range problem 
has oscillatory solutions when a positive margin of 
thrust over drag exists. With zero thrust the 
solution obtained is the flattest glide with a 
family of transients to and from the maximum lift- 
to-drag point. For the time-range problem, values 
of H/Xx to the left (low speed end) of the V(T-D)max 
point produce oscillatory solution while on the 
right of the V(T-D)m=ix point a family of transients 
to and from the equilibrium point defined by the 
choice of H/Xx is obtained. 

Leqendre-Clebsch necessary condition ; 

From the Euler-Lagrange equations, with Xj^ = 0 

3H = -Xw g Cosy + V Cosy - X. V Siny (54) 

3y ^ 

and 


O. = ( Xy g - X^ V ) Siny - V Cosy 


(55) 


3y" 

Setting the left-hand side of equation (54) to 
zero as required for a stationary minimum of H leads 
to 

Tany = X^V - Xyg or Xyg - X^V (56) 






From (56), then Siny = ( X|^V - Xyg )a 


s/CTJ 




(57) 
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and 

Cosy " 




- V 


(56) 


where j ■ + 1 


Using (57) and (58) in (55), it Is possible to 
determine 3. 


first valjje of the Independent variable x » x* > Xq 
for which such a nontrivial system can be found 
defines a conjugate point. 

Following the analysis of Ref. 9 for the Mayer 
problem, the rank of the matrix of variations of 
states and the multiplier corresponding to the state 
being minimized with respect to the Initial values 
of costates Is evaluated along the test extremal, 

VIZ. 


Next, one may employ the transversal Ity 
conditions for the range problem, vlz.x^ » 1 for 
range minimization and » -1 for range maximiza- 
tion. These lead to 

« 1, 3^H >0 Ify lies In the second or third 
quadrant (S®) 

Xjj =-l, ifit <0 If Y lies in the first or fourth 
3 y^ quadrant (60) 

It is apparent that with no restrictions on 
the path Inclination, the minimum-range trajectory 
Is that which maximizes the range In the negative 
direction. If one imposes an artificial limit on 
the path angle y, say -go y ^ + 90, the minimum- 
range trajectories are vertical up-down flight 
segments and any (V,h) pair can be reached In zero 
range If there exists a positive margin of thrust 
over drag. Thus, with the vehicle model considered, 
there Is no limit to the steepness of climb. 

The Legendre-Clebsch necessary conditions are 
met 1n strengthened form for the maximum-range 
problem for values of the path angle y In the first 
or fourth quadrants. However, physical reasoning 
makes clear that a range maximization problem with- 
out time or fuel constraints will not possess a 
proper maximum or even an upper bound. In view of 
the above, the problem of Interest is to maximize 
the range to climb from an Initial (V,h) pair to a 
final (V,h) pair In a fixed time. This problem is 
of value in studies of the type reported in Ref. 8 
for a point-mass-modelled vehicle. 

It may be noted at this point that in the cases 
of time and/or fuel minimization problems with 
range open, the Legendre-Clebsch necessary condition 
is met only in weak form along central arcs and, 
hence, these trajectories fall into the class of 
singular extremals, 

Conjugate-point test : 

The Legendre-Clebsch necessary condition is met 
with a margin for the time-range problem and hence 
the Euler solution (17) with = 0 furnishes a 
relative minimum for initial and terminal points 
sufficiently close together. For extremals of 
finite length, however, the task of ensuring that 
the second variation is non-negative for admissible 
neighboring paths leads to the accessory-minimum 
problem in the calculus of variations. This in 
essence bolls down to a search for a system of ad- 
missible variations, not Identically zero, which 
offer the most severe competition in the sense of 
minimizing the second variation, If a system of 
nonzero variations can be found which makes the 
second variation zero, then it 1s clear that a 
neighboring path •’’s competitive and that the tes-*- 
extremal furnishes at best an improper minimum and 
at worst a merely stationary value (Ref, 9), The 


The rank of 



3 X 2 ^ ^ ^ 3 X 2 



( 61 ) . 


3X„ 

3^n 

* « 1 » * n 

3X, 

Iq 

3X, 


3X, 

3X, 


1 n 

3^ ^ 



provides the criterion for the existence of a conju- 
gate point. If the rank of the test matrix (61) 
drops at any point along the test extremal, it is 
indicative of the occurrence of a conjugate point. 


For the time-range problem, If the Independent 
variable Is changed from time to range, the equa- 
tions of motion become 

h = Tanv (®^) 

v' = alML. -iTanY (63) 

WV Cosy v ' 

The optimal -control problem then is to maximize 
the final value of altitude *h' for a specified 
range with time fixed. With the interpretation of 
H as the time multiplier, the test matrix (61) be- 
comes 


3V 

3V 

3V 


3V 

3V 

• 

3V 


3Hq 




3Hg 


3t 0 

3t 

3t 


3t ° 

3t 

3t 


3Hq 

3^0 


3Xu 

^0 

3Hg 

3Yq 


3X^ 

3Xh 


1 

0 

0 

0 

1 

3Hq 

3Yq 






Note that time appears in this problem as a state- 
like variable with 



(65) 


A prime on the variables denote differentiation with 
respect to the range variable x. 

From (64), the sign of 

3V ,3t - 3V -it (66) ■ 

3Yq 3Hq 3Hg 3 Yq 

evaluated along the Euler solution determines the 
rank of the matrix '''-i). If the sign changes at 
any point on the t\ nge trajectory it is indi- 
cative 0 ^ a conjugatb 


The Euler solution obtained for the time-range 
problem, with altitude dependence of u suppressed, 
may now be tested for conjugate points, In view of 
the particularly simple form of the conjugate-point- 
test for this problem, 1t seems reasonable to 
attempt to obtain analytical approximations for the 
partial derivatives In (66). 

Linearizing the equations of motion and the 
Euler equation (17) with range as the Independent 
variable about an equilibrium point at a particular 
altitude, 

6V' = 3q - a, ay (67) 


T « 


" ^3 


ORfGirv'At PACfi 

OF POOR QUAUTV 


( 81 ) 


and cancelling out coiuion constants In the numera- 
tor, the e ’uatlons (76) - (79) can be brought to 
the form 



n 

JTT 


S^ -toj' 


n 


(82) 



s(s^+.i^) 


(83) 


at * -ag av -t- aj ay 

ay = a^ av - agay ag aH 


(68) , (HTs)..^- 

(69) s(sV„) 


(84) 


Where : 
a 


’0 


q 3 (T-D) 

RVfe W 


'1 


= 1 

^ “5 

rcosy 

a-j = Siny 

0 o 

VCos'^y 

3;, 5 -Cosy q H 3(T-D) 


(70) 

(71) 

(72) 


X 

+ _2 


V^(T-D) 


t 


Cos^ 

V 


^ .llMfri - Cosx h1 

V(T-D)^ \ [ V 




vTt^ 


3^(T-D ) 


3V^ 


Cosy H _ 1 


^5 " 


3^1-D ) + (T-0) j 


(75) 


36 6os_y_.. a. 

V^'T-0) ^x 

Equations (67), (68) and (69) constitute a linear, 
constant-coefficient system which can be put in the 
following form using Laplace transforms (initial 
conditions on and it are zero) 




mm- 7i?^) 


(85) 


Equations (84) and (85) may be further simpli- 
fied using the expression (82) and (83), 

atls) => + T ayM (86) 

ay( 0 ) 


ay 


aiil 


at(s ) = 
W) s 


W 


+ T av(s! 

afiw 


Equations (86) and (87) imply 
at(x ) = §Ul) + T aV(x ) 

X., XU ' Sv- 


ayg 


ayg 


(73) 

at(x) = 

1 

CSJ C 
3 

(74) 

6Hg 



+ T ayix) 
aHg 


(87) 

( 88 ) 
(89) 


using (88) and (89) in (66) 

- 3V . at i av(x)|-i _ ***n 
^0 ^ 


3yg ^ 


■n'J 



(90 


oV(s) = 

■ ®1 

OT) 

s^+(a-|a4-ag35) 

oV(s) = 

■®1^6 

Ho) 

s[s2+(a^a4-aoag)] 

5t(s) = 


aynr) 

s [s^+(a^a 4 -aga 5 )] 

dt(s) = 

■ [(3o^3*®1^2^‘®3^ 3 ®6 

HT) 

s2[s^+(a.,a4«-aQa5)] 

putting 

= (ai>^4"^0®5^ 


and 


(76) 

(77) 

(78) 

(79) 

(80) 


And consequently, one needs to obtain the inverse 
transform of only three transfer functions, namely 



av(s ) . y ^ 

w) 7(7^} 


2 

when % is positive, the roots of the denominator 
polynomial are conjugates and 

aV . n - ^ n * u xSin((o„x) -f 2Cos(u^x)- 2 

3yg 3Hg 3Hg 3yg (91) 


The right hand side of (91), after being zero at 
X = 0, will subsequently become zero at 


X 


= ^ 

CU 

n 


(92) 


implying that conjugate points will occur every full 
cycle of oscillatory solution. Hence, if the 
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equnibHum point for the given H/^x Is stable. I.e. 
It produces an oscillatory solution, a conjugate 
point will occur at the end of one full cycle of 
the oscillation. On the other hand. If w,j* 1$ 
negative, the roots are real and distinct, symmetric 
about the Imaginary axis. In this case 


, n. - 3V at * 

3Y0 aFg 5yg 


(93) 


-x.d.Slnh(dx) + 2 Cosh(dx) - 2 


Where 



Expression (93) Is zero only at x » 0. In this case, 
conjugate points do not occur. From (93), then, If 
the equilibrium po’int for the given H/^x 'Is unstable, 
conjugate points will not occur. 


The conjugate-point test Is now applied to the 
three regimes of H/\x described earlier. As ex- 
pected, for all values of H/X^ to the left of V(T-D) 
max point, conjugate points occur. Indicating that 
the Euler solutions obtained with these values of 
H/'x do not afford a maximum to the time-range 
problem over long Intervals. Euler solutions ob- 
tained with H/xx to the right of the V(T-0)max point, 
on the other hand, satisfy the Legendre-Clebsch 
necessary conditions and Jacobi's necessary con- 
dition, and hence are optimal trajectories for the 
time-range problem. 


Numerical solution of the time-range problem ; 


With the insight gained for the time-range 
problem with altitude dependence of thrust and drag 
suppressed, we now embark upon a numerical study of 
the more general case 1n which the aerodynamic coef- 
ficients are functions of Mach number and the thrust 
is Mach-altitude dependent. The data for a version 
of the F-4 aircraft with afterburner operative are 
used in this study, A cubic-spline representation 
(Ref, 10) Is used to compute the values of zero- 
11 ft drag coefficient and the induced-drag coef- 
ficient, The drag coefficient is then computed as 

C„ . C„^(M) t K(H) cl 

Where = W and Cg and K are standard 
notati8n. 

The drag is then obtained as the usual product 
of drag coefficient, dynamic pressure and the air- 
craft wing area, A cubic-spline lattice (Ref, 10) 
is used to compute the value of thrust at a given 
altitude and Mach number. Atmosphere density and 
speed of sound as functions of altitude are inter- 
polated from standard atmosphere tables using cubic 
splines. The system differential equations are 
integrated using a fifth-order Runge-Kutta-Verner 
method with variable step size. 

A plot of H/Xx vs airspeed for equilibrium 
flight conditions corresponding to unaccelerated 

v2 

flight with specific energy E = h + ^ , frozen 

at 60,000* is shown in Fig. 9. The three regimes of 
H/'x identified earlier in this paper can be seen in 
Fig. 9. Numerical integration of the Euler equation 
with H/\x values picked from each of these regimes 
Indicated that the solution for H/'x values to the 


left of V(T-D)tnax are oscillatory, Numerical 
solution using H/Xx to the right of V(T-D)mav point 
(high speed end) are non oscillatory and violent In 
character, 

Next, a numerical conjugate-point test Is set 
up based on a scheme suggested by Cicala (Ref. 11). 
In this scheme the partial derivatives with respect 
to 'ig f’equif’ed 1n the matrix (64) are calculated 

approximately In terms of difference quotients. 

Small Increments In Initial x^ are employed In the 
ev/ 4 luat 1 on of neighboring solutions of the original 
system of Euler equation. The conjugate-point test 
was carried out for various values of H/Xx picked 
from Fig. 9. It was found as expected that only the 
non-osclllatory trajectories corresponding to H/Xx 

values on the right of V(T-D);r)ax satisfy the no 
conjugate point condition. Oscillatory trajectories 
Indicate the existence of a conjugate point after a 
cycle of oscillation. 

From the foregoing It is clear that the solu- 
tion to time-range optimal -control problems are 
nonoscillatory and violently unstable In character. 
Within the permissible range of H/Xx, as H/Xx In- 
creases, the Euler solutions approach the energy 
climb schedule 1n the V,h plane. Of particular 
Interest In practical applications Is that trajec- 
tory which terminates at the "dash-point" on the 
flight envelope, the maximum-level-flight-speed 
point. To determine the value of H/Xx which will 
accomplish this, a plot of the locus of equilibrium 
points corresponding to unaccelerated flight at 
constant energy Is made. Once this value H/Xx Is 
found, what remains to obtain the optimal trajectory 
Is to determine the Initial value of the control 
variable y for a given set of initial conditions on 
altitude and airspeed. 

In Fig, 10 the level-flight envelope for the 
F-4 aircraft is shown along with the energy-climb 
schedule, The discontinuity in the energy-climb 
schedule due to transonic drag rise may be noted 
(Ref. 12). The curve B Is the locus of equilibrium 
points at each energy level corresponding to un- 
accelerated flight with the appropriate H/Xx. The 
discontinuity due to transonic drag rise Is again 
visible. An Euler solution for initial values of 
airspeed and altitude close to the equilibrium locus 
is also shown. To determine this trajectory an 
Iteration was undertaken on the Initial value of 
the control variable y. With quadruple precision 
on the lBM-370/158, the initial path angle had to 
be determined to 13 significant digits, To Illus- 
trate the sensitivity of the Euler solution to the 
Initial value of path angle y, the last digit of yq 
is perturbed In the positive and negative sense with 
the trajectories 1 and 2 shown in Fig. 10 resulting. 

A few more Euler solutions with initial condi- 
tions far removed from the equilibrium locus are 
shown in Fig. 11 . 


Discussions and Conclusions; 


In this paper, optimal flight in the vertical 
plane with a vehicle model intermediate in com- 
plexity between point-mass and energy models was 
studied. Flight path angle takes on the role of 
control variable in the model and range-open prob- 
lems feature subarcs of vertical flight and singular 
subarcs as previously studied. 
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Mlnlniuiu-range cUmb problem (the steepest 
of Ref. 15) has been found to have no minimum, not 
even a lower bound. In reference 15, the steepest- 
cllmb problem was studied using the Green's theorem 
device of reference 13 and 14, There Is an Impor- 
tant difference In vehicle modelling from that of 
the present paper which should be noted as a key to 
resolving disparities between the character of 
optimal paths emerging. The analysis of Ref, 13 and 
15 In essence replace cosy In equation (3) with 
unity so that the problem solved is maximum altitude 
In a given distance (I.e.arc length) rather than In 
a given range. This Is a necessity with the linear- 
integral approach which can accomodate only problems 
of dimension two and very special form of state 
equations. The solution to the distance climb con- 
sists of a central path flown along a(T-D)niax locus 
In theV-h plane with vertical climb and dive tran- 
sitions at the ends to meet specified boundary con- 
ditions. 

From physical considerations It can be seen that 
when a positive margin of thrust over drag exists, 
the maximum-range climb trajectory without time or 
fuel constraints has no proper maximum nor an upper 
bound, In view of this fact major attention has 
been accorded to the time-range problem. 

For the special case In which the thrust and 
drag depend only on airspeed, a plot of the ratio 
of time and range multipliers H/x^ for equilibrium, 
corresponding to unaccelerated flight, revealed the 
existence of three regimes. Positive values of 
H/'x tbe low-speed side of V(T-D)fnax 2 ll 
negative values of H/X,^ were shown to yield oscilla- 
tory solutions. Although these meet the Legendre- 
Clebsch necessary conditions, they fall the con- 
jugate-point test. Euler solutions with H/Xv chosen 
to the right of the V(T-D)max Pcfnt satisfy both 
Legendre-Clebsch and Jacobi necessary conditions 
and are nonosclllatory in character. Depending on the 
nature of aircraft data, unstable equilibrium points 
may sometimes appear for certain H/xx values to the 
left of airspeed corresponding to ( T-D at cert- 
ain energy levels, These normally have short duration 
and are not of major interest. 

Numerical solution of the Euler equation and 
a numerical conjugate-point test for the F-4 
aircraft data reinforced the conclusions arrived at 
in the analytical exercise. 

From a practical viewpoint, the time-range tra- 
jectories which terminate at the "dash-point" on 
the level flight envelope are of particular interest. 
The multiplier ratio H/X^ corresponding to this 
point is determined using the locus of equilibrium 
points at each energy level corresponding to un- 
accelerated flight. With this value of H/'x, Euler 
solution for any (h,V) pair is obtained by 
iterating on the initial value of v, the control 
variable. 

Euler solutions were obtained for various 
initial conditions, One observes that these tend 
to funnel rapidly into a certain corridor in the 
V-h chart, in the vicinity of the equilibrium locus 
corresponding to unaccelerated flight. This feature 
of the solution family can be exploited in practical 
situations to simplify the computation of optimal 
trajectories, 
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Fig. 2 . Path Angle vs.Airspeed in Gilding Flight 
for Range Problem 





Fig. 4 . H/X^ vs.Aihspeed at Equilibrium points 

( Parabolic {T-0)/W distribution ) 
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Abstract 

MlnlMwn-tlme climbs In "energy" approximations 
are reviewed and further consideration given to 
choice of variables, A pair of variables which 
seem to offer attractive replacements for altitude 
and airspeed In singular-perturbation procedures Is 
suggested. Use of the new variables In an energy- 
modelled climb-dash problem Is Illustrated, 


ft * V 

(4) 


(5) 

where 0(h,V) Is evaluated for L ■ W. 


"Slow"-Vai’lable Choice 



Introduction 


In Ref, 1, Fritz Kaiser, a flight-test engineer 
at Hesserschmltt, A. G., Introduced the concept of 
"GESAMTHOHE" ("resultant height") In connection 
with aircraft mlnlmum-tlme climbs. This Is the sum 
of potential and kinetic energy per unit weight. 
Subsequently It has been referred to as "energy 
helght"‘»3 and "specific energy". 4 Its use as a 
state variable In trajectory work Is attractive be- 
cause It Is a "slower" variable than either alti- 
tude or velocity.®'®*' Attempts to synthesize 
"slow" state variables are described In Refs. '6 and 
7 In connection with singular-perturbation proce- 
dures. The present development attempts to synthe- 
size both "fast" and "slow" variables for the 
mlnlmum-tlme-to-cllmb problem along lines explored 
earlier In an appendix to Ref. 6. 


Climb Equations 

The equations of motion for climbing flight are 
given In terms of conventional state variables, 
altitude, h, flight-path angle, y, and velocity, 

V, as 


fi = V sin Y 

(1) 

Y f [g ' cos y| 

(2) 

V = - g sin Y 

(3) 


Two new variables, a 'nd are to be Intro- 
duced In place of h and V, <1* to be "slow" and 
"fast". 

The equation of state for Is 



. 1^ . [v t . , |^]v 

(6) 

If one Insists that 4> be Independent of the 
control -like variable, y, then ij> must satisfy the 
partial differential equation 

V If - g H » 0 

(7) 

This Is satisfied by 


w2 

.j, “ h + 2^ 

(8) 


or by any once-dlfferentlable function of this ex- 
pression,® Thus <{i = E, specific energy. Is "slow" 
In the sense specified. 


"Fast"-Var1able-Cholce Considerations 


It has been usual to adopt as the second state 
variable, i|/, either V® or h.' Either Is suitable 
for analysis of the "slow" motion, given by the 
single state equation 


Here T Is thrust, D drag, L lift and g the accele- 
ration of gravity. An assumption of thrust-along- 
the-path has been Incorporated. 


Energy -Modelling Simplifications 

An essential feature of "energy" approximation 
is that drag be treated as a function of h and V 
only. This Is consistent with approximation of 
sin Y and cos y via expansion In powers of y 
through, first-order terms only and with deletion 
of the Y term as negligible -- another feature es- 
sential to reduction In order. With these simpli- 
fications the system becomes 


f (9) 

For minimum-time passage to higher energy levels, 
the right member of (9) Is maximized with respect 
to V or h at constant E. The expression on the 
right of (7) Is "specific excess power", Pe, of the 
flight-performance literature (e.g. Ref, 4; and 
simply p later In the present paper. With a more 
general choice of '('(V,h), the maximization of (9) 

Is done with respect to this variable after V and h 
have been replaced by suitable functions of and 
representing the Inverse transformation. The 
resulting valuej of V, h and E are the same, how- 
ever. 
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Singular-Perturbation Attilys'ls 

The choice of HVih) matters, however, in the 
detenni nation of v along the ”slow-motion'' {or 
**outcr**) solution, os y must be such that ^ » 0, in 
the procedure of Ref, 7, With the choice of » h 
as in Ref, 7, the approximation y * 0 is obtained, 
while if •> V is assumed, then 

Y - (10) 

whicji is, to linear approximation In y, the path 
angle for unaccelerated climb, More generally the 
expression 





» 0 


( 11 ) 


is to be solved to obtain the zeroth-order "outer" 
approximation for y. The choice 



( 12 ) 


suggests Itself for compatibility with the outer 
solution, for 


“ H Vy + If 



(13) 


Here 

(,«) 

is "specific excess power," a known function of h 
and V, This choice of Is seen to generate 
zeroth-order y consistent with (4) and (5) along 
the outer solution. 


Choice of "Fast" Variable 


Contours <ji » E « const, and f » const, are 
shown in Fig, 1 for the aircraft data cf Ref. 9 (a 
version of the F-4). The contours of * constant 
indicate a breakdown of one-to-one mapping 
associated with Jumps of the energy-climb path, 
i/) “ 0, between ridges of p(h,V)^; In fact, the 
mapping -► h, V Is two-to-one and even three-to- 
one within the flight envelope. This local non- 
invert 1 bill ty represents a less -than -Ideal feature 
for a coordinate transformation; however, one does 
not actually have to transform to the new variables 
to exploit the concept, 

Flight-path angle y Is shown as a function of 
^ ■ E in Fig. 2 for three choices of "fast" 
variable: h, V and ii>. The contribution of the 
"outer" solution (shown) dominates the zeroth-order 
composite solution. Experience is that the calcu- 
lation of first- and higher-order composites Is 
Intricate and expensive. 10 Thus It makes sense to 
choose variables carefully so as to enhance the 
fidelity of the zeroth-order solution as far as 
possible. 


Consider as an application the climb-dash 
problem, In which a minimum-time trajectory to a 
reniote value of x is sought, wltere x is down-range 
and, for small y Is defined by 

X - V (15) 

The character of the solution is that of a combined 
climb-dash generally faster than an energy climb 
(Fig. 3) fairing Into sustained flight at the high- 
speed point on the 1ev<e1 -flight envelope, y as a 
function of E is shown in Fig, 4 for the three 
choices of fast variable. Solutions of a cor- 
responding point-mass-modelled problem for dif- 
ferent aircraft data are studied In Ref. 8, 
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Fig.2 0 vs. E- Energy Climb 



E(Feet) . XK)3 

Fig.4 9 vs.E- Climb-Dash 
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Fig.3 Altitude vs. Velocity- Climb- Dash 







